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We investigate the emergence of three-dimensional behavior in a reduced-dimension Bose-Einstein 
condensate trapped by a highly anisotropic potential. We handle the problem analytically by per- 
forming a perturbative Schmidt decomposition of the condensate wave function between the tightly 
confined (transverse) direction(s) and the loosely confined (longitudinal) direction(s). The perturba- 
tion theory is valid when the nonlinear scattering energy is small compared to the transverse energy 
scales. Our approach provides a straightforward way, first, to derive corrections to the transverse 
and longitudinal wave functions of the reduced-dimension approximation and, second, to calculate 
the amount of entanglement that arises between the transverse and longitudinal spatial directions. 
Numerical integration of the three-dimensional Gross-Pitaevskii equation for different cigar-shaped 
potentials and experimentally accessible parameters reveals good agreement with our analytical 
model even for relatively high nonlinearities. In particular, we show that even for such stronger 
nonlinearities the entanglement remains remarkably small, which allows the condensate to be well 
described by a product wave function that corresponds to a single Schmidt term. 



I. INTRODUCTION 

Three-dimensional Bose-Einstein condensates (BECs) 
confined in highly anisotropic traps are known to ex- 
hibit lower-dimensional behavior when the number of 
condensed atoms is well below a critical value [Ij . Under 
such conditions, the tightly confined dimension(s) can be 
effectively neglected because the characteristic transverse 
energy scale far exceeds the scattering (interaction) en- 
ergy of the atomic cloud. 

Among the applications that take advantage of this 
lower-dimensional character, quantum interferometry 
protocols can be particularly sensitive to the true three- 
dimensional nature of the condensate. In recent work [2], 
simulations of a nonlinear BEC interferometer in highly 
elongated geometries showed significant deviations from 
the quasi-lD model with increasing strength of the non- 
linear scattering interaction. Similar effects have been 
shown to impact the propagation of solitons in quasi-lD 
attractive BECs [HHS] ; this is a potential source of prob- 
lems in implementations of matter-wave interferometers. 

Here we study effects associated with the emergence of 
three-dimensional behavior in reduced-dimension BECs 
trapped by highly anisotropic potentials. We develop a 
perturbative Schmidt decomposition of the condensate 
wave function between the tightly confined (transverse) 
direction(s) and the loosely confined (longitudinal) di- 
rection(s). The perturbation theory is valid so long as 
the nonlinear scattering energy is small compared to the 
transverse energy scales. In contrast to variational meth- 
ods ^i8j, corrections to the reduced-dimension approxi- 
mation are found without relying on any a priori assump- 
tions about the condensate wave function or the shape 
of the trapping potential. Because the perturbation for- 
malism is tied to the Schmidt decomposition, it automat- 
ically encodes information about the entanglement be- 



tween the spatial and longitudinal directions. The dom- 
inant Schmidt term corresponds to the optimal product- 
state approximation to the condensate wave function [91 ; 
within this dominant term, the perturbation formalism 
provides corrections to the lowest-order transverse and 
longitudinal wave functions of the reduced-dimension ap- 
proximation; the main effect is a reshaping of the BEC in 
the tightly confined direction as the strength of the non- 
linear scattering interaction increases. The next Schmidt 
term describes the lowest-order entanglement between 
the transverse and longitudinal directions; the pertur- 
bation formalism allows us to calculate the form and 
amount of this entanglement. 

We compare the results of our perturbation the- 
ory with numerical integration of the three-dimensional 
Gross-Pitaevskii (GP) equation. We study the case of 
*^Rb condensates trapped by cigar-shaped potentials for 
various atom numbers and experimentally accessible pa- 
rameters. By considering different trapping potentials, 
we also investigate the dependence of the results on the 
inhomogeneity of the longitudinal potential. We find sur- 
prisingly good agreement even for relatively high nonlin- 
earities; in particular, even for such stronger nonlinear- 
ities, the entanglement between transverse and longitu- 
dinal directions remains remarkably small, which allows 
the condensate to be well described by a product wave 
function that corresponds to a single Schmidt term. 

We begin our discussion by briefly reviewing in Sec. ITT] 
the mean-field description of quasi-reduced-dimension 
BECs, followed by the perturbative derivation of the 
Schmidt decomposition of the condensate wave function 
in Sec. IIIII Comparison of our perturbative approxima- 
tion with numerical solutions of the time-independent, 



three-dimensional GP equation is presented in Sec. IV 
for the parameters of ^''Rb condensates trapped by dif- 
ferent cigar-shaped potentials and various atom numbers. 



The spatial entanglement between transverse and longi- 



tudinal directions is particularly analyzed in Sec. IV C 
Final remarks are given in Sec. |Vl 



II. REDUCED-DIMENSION APPROXIMATION 
TO A EEC IN A HIGHLY ANISOTROPIC TRAP 

In the mean-field approximation, one describes a con- 
densate of N atoms at zero temperature by a wave func- 
tion ip (normalized to unity) that is determined by the 
time-independent Gross-Pitacvskii equation 



^J.1p 



2M 



V^ + V + g{N - 1 



^, (2.1) 



where V is the external trapping potential, /j, is the 
chemical potential, and g = Airti^a/M is the scattering 
strength determined by the s-wave scattering length a 
and the atomic mass M. For brevity, we generally use 



giN-l) 



(2.2) 



in the following. 

In the case of highly anisotropic potentials, the atomic 
cloud is loosely trapped by a potential VL(r) in d di- 
mensions, referred to as longitudinal (L) dimensions, as 
opposed to the remaining D — 3 ~ d transverse degrees 
of freedom (T), which are tightly confined in a potential 
Vt{p)- If the scattering interaction is sufficiently small 
compared to the transverse energy scale, one can neglect 
the effect of the nonlinear interaction on the atomic trans- 
verse degrees of freedom and hence approximate the con- 
densate wave function by a product wave function, 



^o(p,r) ==^o(p)0(r) 



(2.3) 



where ^o(p) is the ground-state wave function of the bare 
transverse potential and (/)(r) is the solution of the d- 
dimensional, longitudinal GP equation 



'2M 



V?, 



^L(r-)+OTT|0(r-)|^U(r) 



^J'L(, 



which is found by plugging the product ansatz (2.31 into 



the GP equation (2.1 1 and projecting the result onto the 
subspace spanned by £,o- Here ^il = ^j, — Eo is the longitu- 
dinal part of the chemical potential, Eq is the transverse 
ground-state energy, and 



m = jd^pMpt 



(2.5) 



In this reduced-dimension approximation, the trans- 
verse and longitudinal degrees of freedom are decoupled, 
which should hold as long as the number of atoms in the 
condensate is small compared to an (upper) critical atom 
number iVr, defined as the number at which the scatter- 
ing energy becomes comparable to the transverse kinetic 
energy, i.e.. 



9 



{Nt - 1)11 = 



2M 



d^p\V^o\ 



(2.6) 



where 77 = J d^r\tfj\'^ is a measure of the inverse volume 
occupied by the condensate wave function ip. 

As N approaches Nt, one can no longer neglect the 
effects of the scattering interaction on the condensate 
transverse degrees of freedom; as a result, the product 
ansatz (2.3) is no longer a good approximation to the 



three-dimensional wave function. Such effects are respon- 
sible not only for redefining the transverse and longitu- 
dinal wave functions, but also for entangling the spatial 
directions. We show below that these effects can be read- 
ily calculated in the perturbative regime where TV is small 
compared to Nt by performing a perturbative Schmidt 
decomposition of the condensate wave function. 



HI. PERTURBATIVE SCHMIDT 

DECOMPOSITION OF CONDENSATE WAVE 

FUNCTION 



Our goal is to find an approximate solution to Eq. (2.1 ) 



that correctly accounts for the nonlinear effects on the 
tightly confined directions that are neglected by the 
reduced-dimension approximation. Instead of proposing 
an alternative to the product ansatz (2.3), we look for a 



solution to the GP equation in the form of the Schmidt 
decomposition. 



C30 



'4'iP,r) = }_^CnXn{p)<Pn{r) , 
n=0 



(3.1) 



where {Xn} and {4>n} form orthonormal Schmidt bases 
in the transverse and longitudinal directions and the c„ 's 
are the (nonnegative) Schmidt coefficients (the squares c^ 
are the eigenvalues of the marginal transverse and lon- 
gitudinal density operators). The decomposition (3.1) 



is guaranteed to exist, but the Schmidt basis must be 
determined from the GP equation (2.1 1. We can assume 



that the condensate wave function and the Schmidt basis 
functions are real. 

As the deviations from the reduced-dimension approx- 
imation arise in a regime where the scattering interaction 
can be considered as a perturbation to the single-particle 
transverse Hamiltonian, we can solve for the Schmidt de- 
composition in successive orders of a perturbation theory. 
We begin by writing the GP equation in the form 



/iV' = (-ffT + eHL + eg\ip\'- 



V', 



(3.2) 



where Ht{l) — ^(fi^/2M)V^,^N -I- Vt{l) is the trans- 
verse (longitudinal) single-particle Hamiltonian and e is 
a formal perturbation parameter. Notice that due to the 
asymmetry of the trapping potential, the longitudinal 
Hamiltonian H^ and the nonlinear scattering interaction 
are treated as of the same size, both being order e smaller 
than the transverse Hamiltonian. 

As e goes to zero, we expect the solution of the 



GP equation (2.1) to reduce to the product wave func- 



tion (2.3 1, in which the Schmidt decomposition has only 



one term. As e increases, the Schmidt decomposition ac- 
quires additional terms. We are thus motivated to treat 
the Schmidt decomposition formally as a power-series ex- 
pansion in e: 



v^-E 



e Xn 



(3.3) 



In developing the perturbation theory, we find it conve- 
nient to absorb the Schmidt coefficients into the trans- 
verse Schmidt basis functions, which thus satisfy the or- 
thogonality relation 



(XnlXm) = c\8„ 



(3.4) 



The longitudinal Schmidt basis functions arc orthonor- 
mal, 



(0n|0»i 



(3.5) 



The consequences of the normalization of the overall wave 
function and of the orthogonality relations ( 3.4 ) and ( 3.5 ) 
are spelled out in Appendix \K\ 

We also look for the chemical potential as an expansion 
in powers of e, 



m=0 

and similarly for the Schmidt basis functions as 
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.m , 



-^nm 



(3.6) 

(3.7) 
(3.8) 



We seek a solution to the 3D GP equation (3.2 1 to first 



order in e. Thus we are looking for a solution of the form 



A>i{p,r) =[xoo(p)+exoi(p)][0oo(^) 
+ eXio(p)0io(T-) , 



£001 (^)] 



(3.9) 



which is spatially entangled (or nonseparable) . We rel- 
egate the details of the straightforward, but tedious 
derivation of the perturbative equations to Appendix |A] 
and only present the results here. 

The m = terms in the n = Schmidt term cor- 
respond, as expected, to the idealized description of a 
quasi-d-dimensional BEG summarized in Sec.|n} xoo = £.0 
is the ground-state wave function of the bare transverse 
potential, and this means that /xq = Eq is the ground- 
state energy of the transverse trap; ^q q is determined by 
the reduced-dimension GP equation (2.4), here written 
as 



M10OO = (Hl +gTiT4>m)4>m , 



(3.10) 



with the nonlinear interaction renormalized by the aver- 
age of Xqq over itself, 

VT = (xoolXoo) = / d^'pxUp) ■ (3-11) 



Hereafter, for brevity, we usually represent spatial inte- 
grals in terms of bra-ket inner products. The longitudinal 
GP equation (13. 101) also determines the first correction. 



/ii, to the chemical potential. 

There are four first-order corrections to be calculated. 
The functions xoi and (/)oi are the first corrections within 
the n = Schmidt term, i.e., to the transverse and lon- 
gitudinal wave functions xoo and <^00j whereas xio and 
(/)io describe the lowest-order spatial entanglement. The 
transverse functions xoi and xio are determined by the 
linear differential equations 



{^0 - Ht)-^ = gixlo - Vt)xoo = im - HT)-r^ 



where 



(3.12) 



7?i EE ((/-ool^ilo) = / d'^r(l)^oir) 



(3.13) 



is the average of the probability distribution (/)qq 
itself, and 



Avl^ 



(000 1'/'; 



00/ 



vi>o 



(3.14) 



is the variance of (f>QQ 



Notice that the first-order trans- 
verse corrections, xoi and xiOj are driven by inhomo- 
geneities in the zero-order transverse profile Xqq. We 
can write a solution of Eq. (3.12) in terms of the eigen- 



functions and eigenenergies of the transverse Hamilto- 
nian Ht, i.e., Hr^n = E^^n, 



Xoi _ Xio 
Vl A?7l 



-5EeJ^"'^° 



' E„ 



Mo 



(3.15) 



The longitudinal function (/)oi is determined by the 
equation 



{^ll - Hl- 3g77T0oo)0oi 
Here 



35 Tt0oo-M20oo ■ (3.16) 



Tr = 



00 



(6.1^0' 



3\2 



~^ E„ 



>0 



Mo 



(3.17) 



is a coupling parameter, which is determined solely by the 
properties of the transverse trap and which characterizes 
the strength of the coupling of transverse and longitudi- 
nal directions; the quantity tj^/Tt can be thought of as 
the relevant quantification of the transverse energy scale 
as far as the perturbation theory is concerned. Projecting 
Eq. (3.16) onto 0oo gives an expression for the correction 



to the chemical potential, 

M2 = 2~gf^T{4>oi\4o) - ifivl + Ar?i)TT , (3.18) 



which shows that Eq. (3.16) is a linear integro-differcntial 
equation for (f>f)i. 

The remaining longitudinal function, (/)io, is given by 
a trivial algebraic equation, 



^10 



*oo - VL 



AriL 



yoo 



(3.19) 



The first-order longitudinal corrections are driven by in- where po = \/H/Mujt ■ It is easy to see that 
homogeneities in the zero-order longitudinal profile cfiQQ. 



It is an easy matter to derive from Eqs. (3.15) and 



(3.19) that 



(XoilXoo) ^ (xiolXoo) = _s V ^^"1^0^^ ==. -gTr 



hol^fm) 



Aru 



(3.20) 
(3.21) 



At the order we are working, the first two Schmidt coef- 
ficients are given by Cg — (xoolxoo) = 1 and 



c? = (xiolxio) = .g^Ar^i Y^ 



3 \2 






(3.22) 



The only nonlinear equation we have to solve is the 
(differential) longitudinal GP equation (3.10) for (/)oo, but 
we are faced with solving the linear intcgro-differential 
equation ( |3.16[ ) for ipQi. It is easier and more instructive 
to combine these two equations into a single nonlinear 
differential equation for the longitudinal contribution to 
the n = Schmidt term, 

00 = 000 + £001 + 0{e^) . (3.23) 

To do this, we write Eqs. ( 3.10[ ) and (3.16) in the forms 

= (mi - -H^l)0oo - 9Vt<PI + 3e5?7T0o0oi + O(e^) , 

(3.24) 

(3.25) 



= (^1 - i/L)e0oi - 3e5?7T0o0oi 



Identifying /i^ = /ii + e/i2 as the longitudinal part of 
the chemical potential, we can add these two equations 
(and then set e = 1) to obtain a GP-like equation for 0o, 
accurate to first order in e: 

/iL0o = {Hl + 9Vt4>1 - 3.g^TT0o)0o 

= [Hl + gMN - l)4>l " Sg^TriN ~ if^b^] 0o • 

(3.26) 

In the second form of the right-hand side, we restore the 
N dependence to reveal the intrinsic coupling strengths. 
Relative to a GP equation, this longitudinal equation 
has an additional quintic term, which acts as an effec- 
tive three-body, attractive interaction among the atoms. 
This attractive interaction is mediated by the changes in 
the transverse wave function, as evidenced by the appear- 
ance of the (nonnegative) coupling parameter Ty in the 
coupling strength Sg^Tr- Such a self- focusing interaction 
has also been used to study the propagation of solitons 
in attractive quasi- ID condensates 3 5|. 



The coupling constants in Eq. (3.26) can be calculated 



explicitly for a transverse harmonic potential, which we 
use henceforth. In this case, the transverse ground-state 
wave function is the Gaussian 



Xoo(p) 



(^rpg)WI' 



(3.27) 



Vt 



\V2^Po) 



(3.28) 



Moreover, for a pancake {D = I), we have 

1^ 



iU^I) 



— / r ''IT^ 



(3.29) 



.0, 



n odd. 



For a cigar {D = 2), if we use polar coordinates for the 
transverse eigenfunctions, they take the form Cnrm{p, f), 
with Tir and m being radial and azimuthal quantum 
numbers and with the eigenenergies given by En,.m — 
hujxi'inr + \m\ + 1). Then we find that 



(Cn,-m|Coo) — 2 ""'rixSmO ■ 



(3.30) 



It follows from Eq. (3.17) that T^ is given by 



Tt = < 



-^ ln(8 - 4^3) , D = l (pancake). 



2nu)T 3 



(3.31) 



D — 2 (cigar). 



As a result, the coupling constants for a cigar-shaped trap 
{d — 1) are grjT — 2hwTa and ^g^Tr — G/iWTa^ ln(4/3), 
whereas for a quasi-2D pancake (d ~ 2), we obtain gr]T — 
2V2TThujTPoa and Sg'^Tr = 2ATrhujTpW M^ - 4V^)- 



IV. NUMERICAL SIMULATIONS 

The accuracy of the perturbative Schmidt decomposi- 
tion of the condensate wave function can be checked by 
direct comparison to the numerical solution of the three- 
dimensional GP equation (2.1 ) for various atom numbers. 



A. Trap geometry and numerical integrations 

We restrict our comparisons to the case of highly elon- 
gated (cigar-shaped) condensates of ^^Rb atoms in the 
\F = l,mp = —1) hyperfine state, for which a = 
100.4 ap, with ao being the Bohr radius, trapped by po- 
tentials of the form 



V{p,z)^\{Muj\p''-rkz'i) 



(4.1) 



with g = 2, 4, and 10. Such choice of potentials allows 
us to explore how the hardness of the potentials affects 
the results. A hard-walled longitudinal trap corresponds 
to the limit g — >■ cxj. We set the transverse frequency to 
350 Hz, which gives po — 0.6 pvn. 

For the case of a harmonic longitudinal trap ((7 = 2), 
we set the longitudinal frequency to 3.5 Hz and find that 



Nt — 14 000 atoms ^U\. To compare the simulations for 
the different longitudinal power-law potentials, we choose 
the stiffness parameter k so that Nt has the same value 
for the two other values of q; thus all the traps have 
the same one-dimensional regime of atom numbers. We 
define zq = {h'^/Mky^^'^^^'' as a measure of the bare 
ground-state width in the longitudinal direction (zg sim- 
plifies to the analog of po for a harmonic longitudinal 
trap). With these choices, the aspect ratio of the bare 
traps, Po'- zq, is approximately equal to 1 : 10, 1 : 24, and 
1:57 for g = 2, 4, and 10. 



In addition to integrating the 3D GP equation (2.1), 



we also numerically integrate the quasi-lD GP equa- 
tion (3.10) and the perturbative quintic equation (3.26) 



for the trapping potentials (4.1 ) and different atom num- 



bers [12 . The latter two integrations yield the longitudi- 
nal wave functions 4>oq{z) and (poiz) and also determine 
/ii and /2l. 

Given (J)qq(z), we can determine the remaining Schmidt 
functions: (f>io follows trivially from Eq. (3.19), whereas 
the transverse terms can be found from Eqs. (3.15) 



and (3.30), together with gr]T/hujT — 2a, all of which 
yields 



Xo(p)=eoo(p)-ar;i(7V-l)f;|^ 

Tlr — l '^ 

2"'-n^ ' 



Xio 



ip) = -aArjLiN - I) J2 



(4.2) 
(4.3) 



where the functions Cn,o(p) = e '''/2''oL„^(/92/p2)/^pp 
are the m = (azimuthally symmetric) Laguerre- 
Gaussian eigenfunctions for the two-dimensional har- 
monic potential, with energies En^Q — 2hujTnr- 
From the normalization of Eq. (4.3), we get 



Cl 



= V(Xiolxio) = VLi2(l/4)(Af - l)aAr,L , (4.4) 



where we use the polylogarithm function Lis(z) — 



B. Dominant Schmidt term 

In this subsection, we study the dominant (n = 0) 
Schmidt term, using various quantities to compare the 
predictions of the Schmidt perturbation theory with the 
numerical predictions of the 3D GP equation. We also in- 
clude the predictions of the quasi- ID, reduced-dimension 
approximation to determine how significant the pertur- 
bative Schmidt terms are. 

By integrating out the transverse dimensions from the 
numerical solution of the 3D GP equation, we calculate 
the axial (longitudinal) marginal distribution 



nLiz)= / d'p\^ip,z)\' 



(4.5) 



At the order we are working in perturbation theory, this 
marginal distribution is given, according to Eq. (3.9), by 




FIG. 1. (Color online) Longitudinal marginal distribution (in 
harmonic trap units) for a condensate of 1 000 atoms. The 
discrete points are the results of the integration over the trans- 
verse plane of the 3D GP ground-state solution for different 
trap geometries: circles (blue) signify q = 2, squares (black) 
q = 4, and triangles (red) q = 10. The corresponding solid 
lines represent the distribution given by the n = longitudi- 
nal Schmidt function, |(^o(-2)|'^, whereas the dotted lines show 
the unperturbed distribution |(/)oo(z)p. The marginal distri- 
bution is well described by the Schmidt function for all values 
of q. In the case of a harmonic trap (g = 2), the effect of the 
quintic coupling on |0o(^)|'^ is evident, but becomes less pro- 
nounced for higher q, as the longitudinal distribution becomes 
more homogeneous. 



|0o(-2)P, as the contribution from the n = 1 Schmidt 
term is of higher order. In the quasi-lD approximation, 
this marginal distribution is given by |(/)oo(z)p. 

In Fig. n] we plot nL{z) against the distributions 
|^o(-2)P and |0oo(-z)P for a condensate of 1000 atoms 
and the three different values of q. The marginal dis- 
tribution is very well described by the n — Schmidt 
function 4'o{z) for all the potentials (4.1). The q = 2 
case is particularly interesting, for the effect of the cor- 
rection provided by the quintic coupling in Eq. (3.26) 



proves to be quite noticeable due to the inhomogeneity 
of the harmonic potential. As q increases, however, this 
correction becomes less important, because the axial dis- 
tribution becomes more homogenous as a result of the 
more hard- walled and flat-bottomed potentials. 

The better performance of the Schmidt function (f>Q{z) 
over 000 (-z) is also evident in predictions for the chem- 
ical potential /i as a function of the number of atoms 
in the condensate. According to our perturbative ex- 
pansion, the chemical potential is estimated to be /i = 
huJT + fJ-L = f^T + Ml + A*2 , as opposed to the estimate of 
the quasi-lD approximation, /^id — hujT + P-i, these only 
differ by the longitudinal contribution fi2- The difference 
comes directly from the difference between the quasi-lD 
GP equation (3.10) and the perturbative quintic equa- 



tion (3.26), and it is from integrating these two equations 
that we get p.i and pL- We compare, in Fig. [2J the two 
approximations against the chemical potential given by 
the numerical integration of the 3D GP equation. Devi- 
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FIG. 2. (Color online) Chemical potential in units of hur 
as a function of the number of atoms in the condensate for 
longitudinal potentials ( |4.1[ ). Discrete points represent the 
result obtained from the numerical integration of the 3D GP 
equation, with circles (blue) for g = 2, squares (black) for 
g = 4, and triangles (red) for q = 10. The solid lines give the 
corresponding approximation coming from the perturbative 
quintic equation ( |3.26| l, whereas the dotted lines are the esti- 
mates from the quasi- ID GP equation ( |3.10 l. The correction 
introduced by the perturbation theory becomes quite signifi- 
cant as A'' increases above 1 000 and does indeed lead to a bet- 
ter approximation of the 3D numerical results in comparison 
with the quasi-lD model, although the perturbation theory is 
noticeably failing, even for q — 10, as A'^ approaches 5 000. 



FIG. 4. (Color online) Longitudinal marginal distribution 
(in harmonic trap units) for a condensate of 5 000 atoms. 
The discrete points are the results of the integration over the 
transverse plane of the 3D GP ground-state solution for dif- 
ferent trap geometries: circles (blue) designate q — 2, squares 
(black) g = 4, and triangles (red) g = 10. The correspond- 
ing solid lines represent the distribution given by |(/)o(z)P, 
whereas the dotted lines show the unperturbed distribution 
|<^oo(2)P- For this many atoms, the self-focusing quintic 
term clearly over-corrects the unperturbed distribution, sig- 
naling the breakdown of the perturbation theory. This over- 
correction is especially evident for g = 2. For g = 4 and 
g — 10, the over-correction is not as bad, and the perturba- 
tion theory does a reasonably good job, even for this quite 
large number of atoms. 
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FIG. 3. (Color online) Transverse marginal distribution for 
a condensate of 1 000 atoms in harmonic trap units. The 
discrete points are the results of the integration over the axial 
axis of the 3D GP ground-state solution for different trap 
geometries: circles (blue) denote q — 2, squares (black) g = 4, 
and triangles (red) g — 10. The corresponding solid lines 
represent the distribution |xo(p)P, whereas the black dashed 
line represents the bare (Gaussian) distribution |Coo(p)p. For 
all values of g, the marginal distribution is well described 
by the Schmidt function |xo(p)P, which represents a major 
improvement over the bare distribution |^oo(p)P- 



ations from the quasi- ID model are again well captured 
by the perturbation theory, especially for potentials with 
higher values of q. 

We can also integrate out the longitudinal dimension 
to calculate the transverse (radial) marginal distribution 
nxip) = J dz\'ip{p, z)\'^ and compare it with the approx- 
imate transverse distributions |xo(p)P and |Coo(p)P- As 
shown in Fig. [3) the transverse distribution for a con- 
densate of 1 000 atoms is well described by the Schmidt 
perturbation theory for all the potentials (4.1); the per- 



turbation theory correctly accounts for the spreading of 
the condensate in the radial direction. In contrast to the 
axial profile in Fig.jTJ the correction to the bare radial dis- 
tribution is not affected by the inhomogeneity of the lon- 
gitudinal potentials. Instead, as predicted by Eq. (4.2 1, it 



is set by the length ratio arj^ ] the radial distributions are 
nearly the same because tjl varies only slightly among 
the three values of q. 

The success of our perturbation theory emboldens us 
to push it a bit beyond where it really should work. 
Figure |4] plots the longitudinal marginal distribution for 
5 000 atoms, and Fig.[5]plots the transverse marginal dis- 
tribution for 5 000 atoms. In both figures, we can see the 
breakdown of the perturbation theory, although it per- 
forms surprisingly well, especially for higher values of q, 
given that this atom number is more than a third of the 
way to Nt- 




FIG. 5. (Color online) Transverse marginal distribution for a 
condensate of 5 000 atoms in harmonic trap units. The dis- 
crete points are the results of the integration over the axial 
axis of the 3D GP ground-state solution for different trap ge- 
ometries: circles (blue) denote q = 2, squares (black) q = 4, 
and triangles (red) q — 10. The corresponding solid lines rep- 
resent the distribution |xo(p)P, whereas the black dashed line 
represents the bare (Gaussian) distribution |^oo(p)p. For this 
many atoms the perturbation theory starts to break down, as 
the Schmidt function |xo(/9)P over-corrects the unperturbed 
distribution !Coo(p)P- For all values of q, the marginal distri- 
bution is still reasonably well described by the perturbation 
theory. 



C. Spatial entanglement 

The last two Schmidt functions, Xioip) ^-nd (f>iQ{z), 
introduce the final correction to the quasi- ID, reduced- 
dimension approximation, namely, to the assumption of a 
spatially separable three-dimensional wave function. The 
validity and importance of these corrections can be as- 
sessed in terms of two quantities. The first of these is the 
probability Pjj of finding the numerically determined so- 
lution of the 3D GP equation (3.2 1 outside the subspace 
spanned by the Schmidt wave functions Xo^o and Xi0i- 



This probability deficit tells us the extent to which the ex- 
act solution is confined to the two-dimensional subspace 
of the perturbative Schmidt wave function (3.9 1 and is 
given by 



Pd = 1-?4- 



(4.6) 



where cq and ci are Schmidt-like coefficients obtained by 
projecting the exact 3D solution onto the (normalized) 
Schmidt basis functions, Xo^o/co and xi't'i/ci- Com- 
puted values of the deficit P]j for the different longi- 
tudinal potentials are displayed in Table [l] for various 
atom numbers. The very small values of the deficit in- 
dicate the success of our perturbation theory. For larger 
atom numbers, however, the population outside the two- 
dimensional space increases as the perturbation theory 
begins to break down. 

Within this two-dimensional subspace, we can assess 
the validity of the perturbative wave function in terms of 
an entanglement measure. Since the perturbative wave 



TABLE I. Probability deficit Pd x lO". Small values are in- 
dicative of the overall success of the two-Schmidt-term per- 
turbation theory. 



q 


Af = 1 000 


2 000 


3 000 


4 000 


5 000 


2 


0.19 


1.76 


6.88 


18.05 


39.23 


4 


0.03 


0.44 


1.81 


4.91 


10.41 


10 


0.05 


0.21 


0.88 


1.44 


3.06 



function has the form of a two-qubit entangled state, we 
can use Wootters's concurrence ^13j for a pair of qubits 
as the entanglement measure. 

The concurrence of a bipartite pure state |^as) varies 
smoothly from for product states to 1 for maximally 
entangled states. From its definition, C = K^abIcj/ ^ 



<^yl*Ai3> 



in terms of the Pauli matrix ay and the com- 



plex conjugate of \'^ab), it is easy to show that the 
concurrence for the perturbative condensate wave func- 



tion (3.9 1 is given by 



C - 2coCi = 2v/Li2(l/4)(iV - l)aATjL 



(4.7) 



where we use the perturbative coefficients cq — 1 and ci 
as given by Eq. (4.4). The concurrence tells us about 



the amount of entanglement generated by the nonlin- 
ear interaction between the radial and axial directions. 
This information, which quantifies the importance of the 
nonseparable corrections, is essentially contained in the 
Schmidt coefficient ci. 

In Fig. [61 we compare the concurrence of the exac t 3D 
solution, given by C = 2coCi, with the concurrence (4.7) 
of the perturbative Schmidt theory [14] . One can see 
that the entanglement remains remarkably small even 
for relatively large atom numbers. Notice that as the 
inhomogeneity of the longitudinal potential decreases, so 
does the spatial entanglement, which is an immediate 
consequence of Eq. (3.19). In fact, for homogeneous lon- 



gitudinal potentials, such as rings or boxes, the spatial 
entanglement vanishes completely. 

The probability (4.6) and concurrence play comple- 



mentary roles in assessing the accuracy of the pertur- 
bation theory: the deficit Pjj tells us to what extent the 
exact solution is confined to the two-dimensional sub- 
space of the perturbation theory, and the concurrence C 
tells us whether the perturbation theory captures cor- 
rectly the entanglement in this subspace. 

We can also use another, more directly physical quan- 
tity to quantify the effects of the nonseparable correction, 
namely, the condensate average density Nr] = N{^p\^^). 
This parameter is of special interest in interferometric 
applications of BECs f2]. The average density is par- 
ticularly appealing, because it can be used to measure 
any of the perturbative corrections introduced by the 
Schmidt decomposition, not just the nonseparable cor- 
rections. In Fig. [Tj we plot the average density obtained 
from the three-dimensional GP ground-state solution, as 
well as its estimates obtained from (i) the entire per- 
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FIG. 6. (Color online) Concurrence of the condensate wave 
function as a function of the number of atoms in the con- 
densate. The discrete points are the concurrence calculated 
from the Schmidt coefficients found from the numerical so- 
lution of the 3D GP equation, as described in the text: cir- 
cles (blue) signify q = 2, squares (black) g = 4, and trian- 
gles (red) q = 10. The corresponding lines are the concur- 
rence of the perturbative Schmidt wave function, as given by 
Eq. ( |4.7[ ). The entanglement between transverse and longi- 
tudinal dimensions decreases for higher values of q, since the 
potential becomes more homogeneous. Moreover, the entan- 
glement remains remarkably small even for relatively large 
A'^, indicating that the condensate is well approximated by a 
single Schmidt term. 



turbative Schmidt wave function (3.9), (11) the domi- 
nant Schmidt term alone, Xo^Oj and (ill) the quasi-lD, 
reduced-dimension approximation, which gives average 
density NrjLrjT- Notice that the approximation provided 
by the Schmidt decomposition is remarkably good even 
for relatively large atom numbers. The same is true for 
the approximation given by the dominant Schmidt term, 
which demonstrates how small is the effect of the nonsep- 
arable corrections. Not surprisingly, in view of our pre- 
vious results, the unperturbed estimate Nrj^riT quickly 
deviates from the exact numerical results. 

The results of this subsection point to the interesting 
conclusion that the nonseparable corrections to the dom- 
inant Schmidt term are small even for relatively large N . 
Therefore, as far as the ground-state properties go, one 
can nearly neglect the spatial entanglement and describe 
the condensate wave function in terms of the product 
wave function corresponding to the first Schmidt term. 



V. CONCLUSION 

We begin our concluding section by noting that there is 
another way to do the perturbation theory, which we call 
the relative-state method, in contrast to the Schmidt- 
decomposition method developed in this paper. The 
relative-state method starts from a relative-state decom- 
position of the condensate wave function in the bare 
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FIG. 7. (Color online) Average atomic density N-q in har- 
monic trap units as a function of the number of atoms in the 
condensate. Discrete points are the results obtained from the 
numerical solution of the 3D GP equation for different trap 
geometries: circles (blue) are for g = 2, squares (black) for 
q = 4, and triangles (red) for q — 10. The corresponding solid 
lines are those obtained from the entire perturbative Schmidt 
wave function (3.91, and the dashed lines are those of the 



dominant Schmidt term only, whereas the dotted lines, which 
deviate substantially from the exact numerical results, are 
the result, NtitVl, obtained from the quasi-lD approxima- 
tion. The agreement between the exact results and those of 
the perturbation theory, even for relatively large atom num- 
bers, is remarkable. The difference between the solid and 
dashed lines comes from the nonseparable corrections to the 
dominant Schmidt term; this contribution is quite small, con- 
firming the conclusions drawn from the concurrence. 



transverse eigenbasis ^„: 

oo 

V'(P, r) = ^a{p)ipo{r) + e ^ ^nip)^n{r) . (5.1) 

ri=l 

The dominant longitudinal wave function, (/jq, is de- 



termined by projecting the 3D GP equation (2.1) onto 



fo and thus satisfies the reduced-dimension GP equa- 
tion (2.4); it is not the same as (pQ of the Schmidt the- 
ory, but rather is the same as ^oo- The relative-state 
longitudinal wave functions (/?„ are neither orthogonal 
nor normalized. One can work out the relative-state 
perturbation theory order by order, and this is more 
straightforward than the same process in the Schmidt- 
decomposition perturbation theory. At the order we are 
working, we have shown explicitly that the two pertur- 
bation theories are equivalent. 

Even though the relative-state perturbation theory is 
easier to implement, especially at higher orders than we 
consider in this paper, it obscures the physics of the con- 
densate wave function. In the Schmidt perturbation the- 
ory, the n = term gives the best product approxima- 
tion to the exact wave function, and the Schmidt theory 
gets directly at how that term is modified by higher- 
order effects of the nonlinear scattering interaction. The 
n > terms in the Schmidt theory describe entangle- 
ment of the transverse and longitudinal directions, and 
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the Schmidt approach neatly separates this entanglement 
from the dominant (n = 0) Schmidt term. The relative- 
state theory can be used to calculate all these effects — 
and, as noted, it gives the same results at the order we 
are working — but it does not divide up the terms in per- 
turbation theory in this neatly interpreted way. 

We close with two ideas, on both of which we have al- 
ready commenced work. The comparisons of the Schmidt 
perturbation theory with the exact numerical results, 
presented in Sec. |IV[ suggest that the ground-state con- 
densate wave function, for surprisingly high atom num- 
bers, is well approximated by the two-Schmidt-term per- 
turbation theory and even by the dominant Schmidt 
term alone. This prompts us to wonder what happens 
to the Schmidt decomposition of the ground-state wave 
function as the atom number crosses over from reduced- 
dimension behavior for atom numbers below Nt to the 
three-dimensional behavior above Nt- We have evi- 
dence, to be presented elsewhere, that throughout the 
cross-over, the ground-state wave function continues to 
be well approximated by two Schmidt terms, and this 
suggests that an approximate description based on just 
two Schmidt terms might work through the entire cross- 
over and into the three-dimensional regime. 

As we noted out above, the corrections to the reduced- 
dimension approximation are particularly important for 
interferometric schemes that use BECs [2]. A remaining 
problem in that regard is to assess how such corrections 
propagate in time and affect the condensate dynamics 
in a highly anisotropic trap. We plan to address these 
questions in a future publication, which will develop the 



Schmidt perturbation theory for the time-dependent GP 
equation. 
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Appendix A: Perturbation expansion 

Our objective in this Appendix is to derive the equa- 
tions for xoo, 000, Xoi, 001, Xio, and ^lo- The first two 
of these give the lowest-order contribution to the dom- 
inant (n = 0) Schmidt term, and the second two are 
the first-order corrections to the lowest-order behavior 
of this dominant Schmidt term. The last two give the 
lowest-order contribution to the second Schmidt term 
and thus describe the lowest-order entanglement between 
the transverse and longitudinal directions. 

We only need to work to at most second order in e to 
determine the quantities we arc interested in, so we use 
the Schmidt decomposition (3.31 in the form 



V' = XO0O + exi0i + e^X202 + O(e^) 



(Al) 



The transverse and longitudinal basis functions are fur- 
ther expanded as in Eqs. (3.7) and (3.8). 



We first consider the consequences of the various normalization and orthogonality conditions. The overall normal- 
ization of the wave function, 



l = (^IV') = (XolXo)+e'(XilXi)+0(e4) 

= (Xoolxoo) + 2e(xQi|xoo> + e^(2(xo2|xoo> + (xoilxoi) + (xiolxio)) + 0{t^) , 



(A2) 



implies that 



(Xoolxoo) = 1 , 
(Xoilxoo) = , 



(X02|X00> 



(Xoilxoi) + (Xiolxio) 



(A3) 
(A4) 

(A5) 



The orthogonality relations (3.4) for the transverse basis functions, carried to second order in e, imply 

CnSnm = (XnlXm) = (XnolXmo) + ^((XnllXmo) + (X«olXml)) + e^((Xn2|XmO> + (X«olXm2> + (XnllXml)) + O(e^) 



(A6) 



For the cases of interest, we get 



n = m = 
n = l,m = 
n = 2,jn = 

n — m — 1 



Co = (Xolxo) = 1 - e'(xiolxio) + O(e^) , 
(Xiolxoo) == , (xiilxoo) = -(Xiolxoi) , 
(X2o|xoo) =0 , 
c? = (Xilxi) = (Xiolxio) +0(e) . 



(A7) 

(A8) 

(A9) 

(AlO) 
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Likewise, the orthonormality (3.5) of the longitudinal basis functions, carried to the same order in e, gives 



n = m = 
n = l,m — 
n = 2,m ~ 

n = m — 1 



('/'ool'/'oo) = 1 , (0oi|0oo) = , (MM = -{<Poi\(t>oi)/'2 

{4'io\4>oo) = , {4>ii\4'oo) = -'{4'iq\4'oi) , 

{4>2o\4>Oo) = , 
('/'lol^io) = 1 • 



(All) 
(A12) 
(A13) 
(AM) 



Now we use the Schmidt decomposition ( Al I to expand the GP equation (3.2) to second order in powers of e, i.e.. 



(^o+eMi + e^M2)xo0o + (a^o + eMi)eXi0i + Moe^X202 

= Ht{xo<J>o + eXi-^i + e'x202) + ({Hl + .9Xo'/'o)Xo0o + e'l^L + 3.gXo'/'o)Xi'^i + O(e^) 



(A15) 



basis functions, we get 



By projecting Eq. ( A15 ) onto (f>o and then onto xo, keeping in mind the strict orthogonality and reality of the Schmidt 



{Ht - Aio)xo + e[(0o|(i?L - Mi)l0o) + 5(0o|'/>o)Xo]xo 

+ e' [ - M2X0 + (MiHL - /ii)|0i)xi + 3ff((/.i|0g)xoXi] + ^(e^) = , 
{Xo\{Ht -Aio)|xo)'?!'o + e[(xo|(-ffT-Aio)|xi)'/'i + (xo|xo)(^l -A*i)'/'o + 9{xo\Xq)<I>o] 
+ e' [ - (xolxo)/i2'/'o + {xo\{Ht - /io)|x2)02 + 35(xi 1X0)00-^1] + 0(6^) = . 



(A16) 



(A17) 



We now expand the Schmidt-basis functions in powers of e as in Eqs. (3.7) and (3.8). Keeping terms to second 



order, we find that the transverse part of the first Schmidt term is determined by 



(A18) 
(A19) 



(^0 - Ht)xqo = , 

ifiQ- Ht)xoi = [{4>qo\{Hl - fJ.i)\(f>oo) +5'7LXoo]xoo , 

(po- Ht)xo2 = [2(001 1 (-ffL -/ii)|(?!'oo) +4g((/)oi|0oo)Xoo~Ai2]xoo 

+ [((PooKHl - Mi)l0oo) + 3otlXoo]xoi + [(0io|(i?L - Mi)|0oo) + 3.g(0io|0oo)Xoo]xio , (A20) 

and the longitudinal part by 

(xoo|(-ffT-Mo)lxoo)'/'oo = 0, (A21) 

[(Xoo|xoo)(a*i - Hl) -5'7t0oo]'/'oo = 2(xoi|(-ffT - Mo)lxoo>0oo + {xoo\{Ht - f^o)\xoo)<Poi 

+ (xio|(i/T-A^o)|xoo)0io, (A22) 

[(xoo|xoo)(/^i - Hl) - 3ott</'oo] '/'oi ^ [{xailiHr - ^J.o)\xol) + 2(xo2|(-ffT - Mo)lxoo) + 2(xoi|xoo)(-ffL - A^i) 

+ 43(xoi|Xoo)<?^oo - (Xoo 1x00)^*2] 000 
+ 2(xoi|(ffT -Mo)lxoo)0oi 

+ [(xoiK^T - /^o)lxio) + (xiiK^T - A*o)|xoo) + 35(xio|Xoo>'^oo]</'io 
+ (xooK-ffr - Aio)|xoo)0o2 + (xioK^t - Mo)|xoq)0ii 
+ (X2o|(i?T - Aio)|xoo)02o , (A23) 



where rjT and jjl-, defined in Eqs. (3.11) and (3.13), are the lowest-order terms in the expansions of (xolXo) ^iid 
(0ol0o)- We now use the normalization and orthogonality conditions from above and the transverse equations to 
discard the first longitudinal equation and to simplify considerably the other two: 



(Mi 



Hr 



9Vt (I'm) 4'oo = , 



(/ii - Hl - 3ott<?!'oo)'/'oi = [(xoil(^T-Mo)lxoi) + 4.9(xoi|Xoo)</'oo -A*2]<?!'oo 

+ [(xoil(^T-Mo)lxio) +3g(xio|Xoo)<^oo]'/'io ■ 



(A24) 



(A25) 



The lowest-order [m — 0) equations for the domina nt {n = 0) Schmidt term are the transverse Schrodinger 
equation (AI8) for xoo ^-^id a longitudinal GP equation (A24) for ^qq. Both xoo ^-nd ^oo are normalized to unity. 



These lowest-order equations are precisely those that give the reduced-dimension approximation discussed in Sec. ^U 
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We can now use the longitudinal equations to simplify the transverse equations and these simplified transverse 
equations, in turn, to simplify further the longitudinal equations, the results being 



(Mo - Ht)xqo = , 

(^0 - Ht)xoi = .m(Xoo - VT)xoa 



(A26) 

(A27) 



{fio - Ht)xo2 = [2g{(j)oi\(t)lo}{2xlo - Vt) - m]xoQ + gVLi^xlo - '7t)xoi + 5('^io|'/'oo)(3Xoo " Vt)xio 



{fii - Hl- 3gr]T(l>lo)(l>oi = [5(xoi|Xoo)(4'/'oo - Vl) - Ai2]'/'oo + g(xio|Xoo)(3<?!'oo " Vl)(I>io 



(A28) 
(A29) 
(A30) 



Notice that the right-hand side of Eq. (A27) is orthogonal to xooi as required by the left-hand side. 

The remaining first-order terms, xw and ^lo, provide the first correction to a product wave function and thus 
describe to lowest order the entanglement of the transverse and longitudinal dimensions. These terms are determined 
by projecting Eq. ( A15) onto (f>i and Xi- The first of these gives 



(Hj 



fJ-oJXi 



(Hr 



Mi)l0o)xo+g(</'i|'/'o)Xo + O(e) 







(A31) 



Plugging in the expansions (3.7) and J3.8L we find that 



{^iQ- Ht)xio = [{(I)io\{Hl - fJ-iMoo) +ff(0io|0oo>Xoo]xoo = ff('/'io|0oo)(Xo - ^t)xoo 



(A32) 



where to obtain the second form, we use the longitudinal equation (A29l for 000- The right-hand side of Eq. (A32) 
is orthogonal to x qq, as required by the left-hand side. 
Projecting Eq. ( |A15| ) onto xi yields 

{xi\{Ht ~ Aio)lxo>0o + e[(xi|(^T - A^o)lxi)0i + ff(Xi|Xo)0o] + 0{e^) = . (A33) 

Applying the expansions (3.7) and ( |3.8| gives at lowest order, (xioK-^^T — Mo)lxoo)</'oo = 0, which is already satisfied, 
and at the next order. 



iXiolif^o - HT)\xio)(f>io = [{Xii\{Ht - fJ-o)\Xoo) + (XioK-^T - Mo)IXoi) + .9(XiolXoo)'/'oo]<^oo 

+ {xio\{Ht - A*o)|xoo)0oi • 



(A34) 



Simpl ifying this equation using Eqs. ( |A26 ), (A27), and 
( A32 ) gives a remarkably simple algebraic equation for 



^10, 



»io| 



4) 



^10 



= -^00 - VL(, 



'00 



(A35) 



This means that 0qq is a linear combination of the or- 
thonormal functions (boo and d>To and hence that 



('^ool'/'oo) =vl + ('/'io|0oo> 



(A36) 



•oo|0oo) is the average value of ^qo o'^^'' the 
)qo and {(t>oo\^oo) ^^ ^^e second moment of 
00' ('/'lol^oo)^ is the (nonegative) variance of 0qq, which 

2 



Since rjL = ( 
distribution 



we denote as A?7£. Putting all this together, we have 



^lol^oo) = ^Vl 



(•^ool'^oo)-'?! 



^10 



*go 



^m 



(A37) 
(A38) 



We now let {£,n} be the set of energy eigenfunctions 
of the bare transverse potential, with _E„ denoting the 
corresponding energy eigenvalues, i.e., Hx^n = En^n- 



From Eq. (|A26 ) , we h ave xoo = Co and /iq = Eq- Equa- 
tions (\A27l and (|A32|) imply that 



Xoi 

VL 



Xio 






(A39) 



and this, in turn, gives us the relations (3.20) 



The one remaining step we need to take is to use the 
relations (3.20) and our result (A38l for 0io to simplify 



the equation, Eq. (|A30|), for 
(Ml - Hl 



»0i: 
35'7T</'oo)'/'oi = 35^Tt0oo - M20oo , (A40) 



where T^ is the transverse coupling parameter defined 
by Eq. (3.17). Projecting this equation onto 0oo — or, 



equivalently, projecting the equation for xo2, Eq. (A28l, 



onto Xoo — gives the expression (3.18) for the chemical 



potential. Plugging this expression for /i2 back into the 
equation for (j)oi, Eq. (A40), shows that it is a linear 
intcgro-differential equation for 0oi- Solving this integro- 



differential equation and using the result in Eq. (3.18) 
determines fi2- 

This completes the set of equations we need. Sec- 



tion III summarizes the final forms of our perturbative 
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equations and also derives a GP-like, but quintic equa- 



tion for 00 = 000 + £001: which we use in preference to 
the separate equations for 0oo and 4>oi- 
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